Bulk and wetting phenomena in a colloidal mixture of hard spheres and platelets 
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Density functional theory is used to study binary colloidal fluids consisting of hard spheres and 
thin platelets in their bulk and near a planar hard wall. This system exhibits liquid-liquid coexistence 
of a phase that is rich in spheres (poor in platelets) and a phase that is poor in spheres (rich in 
platelets). For the mixture near a planar hard wall, we find that the phase rich in spheres wets the 
wall completely upon approaching the liquid demixing binodal from the sphere-poor phase, provided 
the concentration of the platelets is smaller than a threshold value which marks a first-order wetting 
transition at coexistence. No layering transitions are found in contrast to recent studies on binary 
mixtures of spheres and non-adsorbing polymers or thin hard rods. 

PACS numbers: 64.70.Ja, 68.08.Bc, 82.70.Dd 



I. INTRODUCTION 



Rich bulk phase diagrams involving colloidal gas, liq- 
uid, and solid phases are found when non-adsorbing poly- 
mers or hard rodlike colloids are added as depletion 
agents to suspensions of colloidal spheres Q, Q- The 
chemical potential of the polymers or rods, with which 
one can tune their concentration, plays a role equiva- 
lent to that of the inverse temperature for a simple one- 
component substance characterized by a soft pair poten- 
tial. Moreover, it has been shown theoretically 0, Q 
and by computer simulation |f| that entropic depletion 
mechanisms lead to interesting wetting phenomena and 
to a rich surface phase behavior in colloidal mixtures of 
spheres and non-interacting polymers or rods. A wetting 
transition and layering transitions have been found when 
the mixtures are exposed to a hard wall. 

Recently the depletion potential between two hard 
spheres due to the presence of hard disclike colloids 
has been investigated d, Q, H> @ ■ Subsequently it has 
been shown within of a free- volume theory that depletion 
induced phase separation in a colloidal sphere-platelet 
mixture should occur at low platelet concentrations in 
systems now experimentally available |lOj |. In view of 
the importance of such suspensions in biomedicine [Tl| 
and geophysics ^3 we investigate in the present pa- 
per bulk and wetting phenomena of sphere-platelet mix- 
tures using density functional theory. We demonstrate 
that the geometry-based density functional theory devel- 
oped for binary mixture of hard spheres and thin rods 
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can be consistently extended to the more 
challenging problem of hard spheres mixed with thin hard 
platelets (Sec. II). Moreover, we study the bulk phase dia- 
gram (Sec. Ill) and the wetting of the mixture at a planar 
hard wall by considering the platelets as thin and non- 
interacting regarding their mutual interactions (Sec. IV). 
Our study provides a direct comparison of the bulk and 
wetting properties of binary sphere-platelet and sphere- 
rod mixtures. 



II. DENSITY FUNCTIONAL AND 
FUNDAMENTAL MEASURE THEORY 

We consider a binary mixture of hard spheres and thin 
circular platelets of radius R s and R p , respectively. The 
number density of the centers of mass of the platelets at 
a point r with an orientation lo p = (9 p ,<j)p) of the nor- 
mal of the platelets is denoted by p p (r,u> p ) while /3 s ( r ) is 
the center-of-mass number density of the spheres. The 
equilibrium density profiles of the inhomogeneous mix- 
ture under the influence of external potentials V e2; t iS (r) 
and V ex t. p (v,uj p ) minimize the grand potential functional 

= J d 3 r p.(r) [k B T (\n[Klp s {v)] - l) - p s + F esM (r)] 
+ ^ J d 3 rduj p p p (r,cj p ) [k B T (ln[A|/> p (r, u) p )] - l) 

- H p + Vext,p(*,Up)] + F ex [p s ,p p ] , (1) 

where A s , A p are the thermal de Broglie wavelengths 
and p s , p p are the chemical potentials of the spheres and 
platelets, respectively. The spatial integrals run over the 
volume V which is accessible to the centers of the par- 
ticles and J duj = J* d6 d<fi. The excess free energy 
functional is obtained by integrating over an excess free 
energy density, 

F ex [p SlPp } = Ml y d 3 rdcu p m4 s \n^,4P\n^}), 

(2) 

where v = 0, 1, 2, 3, i = 1, 2, and r = 0, 1, 2. In Eq. © 
the spatial and angular arguments of the weighted den- 
sities ni s \ n- s \ nr'\ and are suppressed in the no- 
tation. Here we use the following decomposition of the 
excess free energy density $: 

$ = $ s + $ sp , (3) 
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is the original Rosenfeld excess free energy density 
|17j for a pure hard sphere fluid. (For the subtle is- 
sue of the range of validity of the Rosenfeld functional 
at high densities see Refs. 0, &s P takes sphere- 

platelet interactions into account up to first order in the 
number density of the platelets (see discussion below). 
There no is contribution $ pp to the excess free energy $ 
in Eq. (0J because the platelets are considered as non- 
interacting particles regarding their mutual interactions. 
In the present application of density functional theory we 
concentrate on ordering effects induced by a planar hard 
wall such that the resulting density profile of the spheres 
depends on a single spatial variable z in the direction nor- 
mal to the wall. Hence p s (r) — p s {z) apart from possible 
surface freezing at high densities. Moreover, we assume 
invariance with respect to rotations around the z axis by 
an angle <p p , so that p p (r, oj p ) = p p (z, 9 p ), where 9 P is the 
angle between the normal of a platelet and the z-axis (see 
Fig. . In this planar geometry the weighted densities 
are given by 
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p s (z)-kw\ s \z), 

n^(z,e p ) = Pp {z,e p )*wP(z,e p ) 
= Ps (z)*w( sp \z,e p ), 
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where the asterics * denotes the spatial convolution: 
g(z)-kh(z) — J dz\ g(z\)h(z — z\) = 5* h. Note that ri^\ 

(s) (v) 

n t , and rvf are weighted densities which involve only 
variables of either species, while n (sp) is a convolution of 
the sphere density with an orientation-dependent weight 
function, combining characteristics of both species. The 
weight functions of the Rosenfeld excess free energy den- 
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where e z is the unit vector pointing along the z axis and 
0(z) is the Heaviside step function. The Mayer func- 
tion / ss (z) of the interaction potential between two hard 
spheres is obtained through 
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The Mayer function equals —1 if the spheres intersect 
or touch each other and is zero otherwise. The remain- 
ing weight functions can be expressed as (see discussion 
below) 



w[ p \z,9 p ) = 
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These weight functions allow one to generate the Mayer 
function f sp (z) of the interaction potential between a 
hard sphere and a thin platelet through 

-f sp (z,9 p ) = w { p Kwi s) +w[ p) *w { 2 sp) +w[ s) *wi p) . (18) 

Equations - (|18fl completely specify the density func- 
tional theory for the system under consideration. 

Before studying the binary mixture of spheres and thin 
platelets in the bulk and near a hard wall it is instructive 
to compare the fundamental measure theory with the one 
which has been developed rece ntly for a mixture of hard 
spheres and thin rods [TH , ITil ITaL fl6| . For thin rods of 
length L the weight functions corresponding to Eqs. i|15|) 
- <(T7J) are given by 
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FIG. 1: Geometries relevant for the determination of the 
weight functions [Eqs. 1151 . 1161 . and 11911 - 1211 1 for thin cir- 
cular platelets of radius R p [(a) , (b)] and thin rods of length L 
[(c), (d)]. The angle between the normal of a platelet and the 
z-axis is denoted by 9 P while the angle 8 r characterizes the 
orientation of a rod with respect to the z-axis. Only the pro- 
jections of the platelets and rods on the planes of the figures 
are shown. 
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Here 9 r is the angle between the rod and the z axis (see 
Fig.nj. The Mayer function f sr (z) of the interaction po- 
tential between a hard sphere and a thin rod is generated 
through 



fsr {z, 



(r) (a) 



(r) (sr) 
■ W - * IB, 
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The weight functions arc linked with a geometrical repre- 
sentation of the particles which is given in terms of fun- 



damental measures defined as Q = / dz w£ , where j — 
s,p, r labels the species, and A = 0, 1, 2, 3 corresponds to 
the Euler characteristic, integral mean curvature, surface, 
and volume [2(j of the particles. For spheres = 1, 

(;[ s) = R si = AttR 2 , and = AttR^/3, whereas 

(v) 

for thin platelets the volume is very small and Q = 1, 

^ — nRp/A, and Q = ZttR 2 . In the case of thin rods 
both the volume and the surface area are very small and 

(r) (r) 

Q ' = 1, Q ' = L/4. For comparison we note that the in- 
tegral mean curvature of a particle can be obtained from 
the general relation £i = f da (1/Ri + 1 / R-2) / (87r), where 
i?i and i?2 are the principle radii of curvature at the 
point a on the surface and da is a surface element. The 



evaluation of this integral is trivial for thin rods and is 
documented for thin platelets in the appendix of Ref. [2l| . 

From their side view, platelets may be regarded as rods 
(see Figs.nj(a) and (c)) and from their top view (in the di- 
rection of the normal to face) as two-dimensional spheres. 
Therefore the functional forms of the weight functions 
w ( p) (z), w { r) (z) [Eqs. and and w[ p) (z), '(z) 
[Eqs. 1(16(1 and 12U|l ] are similar, while the weight func- 
tion w^f\z) [Eq. 116(1 ] takes into account the surface of 
the platelets (see Fig.[T](b)). 

In order to express the sphere-platelet and sphere-rod 
Mayer functions [Eqs. H18|) and l|23|l ] in terms of spatial 
convolution decompositions, additional weight functions 
«/ sp ) (z,9 p ) and w^ Br '(z, 6 r ) have to be introduced. These 
weight functions contain informations about both species 
of the binary mixtures. Since the explicit expression for 
the additional weight function [Eq. (|22|l ] for the sphere- 
rod mixture is amply discussed in Ref. 0] we do not 
document the derivation of the expression for the corre- 
sponding weig ht function w {sp \z,9 p ) [Eq. (JT7|l] for the 
sphere-platelet mixture but analyze the cases 8 P = and 
9 r = 7r/2 in more detail. In these limits the weight func- 
tions reduce to w^(z,9 p = 0) = w^{z,9 r = tt/2) = 
&\/R 2 — z 2 0(i? s — \ z\). Figure|21displays schematic illus- 
trations of the support of the sphere-platelet and sphere- 
rod Mayer functions which can be calculated analytically 
from Eqs. and lf2Hl : 



f sp (z,6 p = Q) = A sp (z)G(R s - \z\). 



A sp (z) = Tr(y/R 2 -z 2 + R 



and 
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(26) 



A sr (z) = n{R 2 s - z 2 ) + 2L V / R 2 - z 2 .(27) 

Due to the steric interaction, the center of mass of a 
platelet [rod] is excluded from an area in the x — y plane 
of size A sp (zi2) [A sr (zi2)] surrounding a sphere, where 
212 is the distance along the z axis between the center of 
mass of the platelet [rod] and the sphere. Moreover it is 
apparent from the figure that the radius R 2 — z 2 and 
hence the weight function w^ sp ' (z , 0) [w^ sr \z, 7r/2)] are 
characteristic of how a sphere looks from the viewpoint 
of a platelet [rod] . 



III. BULK PHASE DIAGRAM 

Based on the density functional given by Eqs. (JJJ - 1(18(1 
and as a prerequisite for our wall-liquid interface study we 
first study the homogeneous bulk fluid with V ex t, s (^) = 
and V ex t lP (v,uj p ) = in a macroscopic volume V. In 
this case the equilibrium density profiles are constant 
[p s (r) = p s and p p {r,uj p ) = p p ] and the Euler-Lagrange 
equations resulting from the stationary conditions are 
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FIG. 2: Illustrations of the steric interactions of hard spheres 
with thin hard platelets and rods, (a) [(c)] Schematic side 
view of a sphere of radius R s and a platelet of radius R p [rod 
of length L\. The normal of the platelet is parallel to the 
2-axis while the rod is oriented parallel to the x-axis. Only 
the projections of the particles on the planes of the figures are 
shown, (b) [(d)] Due to the steric interaction, the center of 
mass of the platelet [rod] is excluded from the hatched region 
surrounding the sphere in the x — y plane. Here the planes 
of the figures are located at z = 212. From these figures the 
area A sp (zi2) = ir(h + R p ) 2 [A sr {z\2) = irh 2 + 2Lh] with 
h = ^R 2 S - z\ 2 can be inferred (see Eq. J2U [p7|l]1. 



given by dtt[p s , p p ]/dp s = and dfl[p s , p p ]/dp p = 0. In 
the present study we restrict our attention to platelet 
number densities p p R p < 0.2 for which the pure platelet 
fluid and also the mixture are in the isotropic phase and 
hence p s is independent of the orientation of the platelets. 
For comparison, the isotropic-nematic phase transition 
of the pure platelet fluid is first order with coexistence 
densities p P iR p = 0.46 and p p nR p = 0.5 according to a 
computer simulation |2l| . The excess free energy density 
per volume can be expressed as 



k B TV 



P P hi a , 



(28) 



with 




l R 2 R p + 2irR 2 p R s 



2(1 -Vs) 2 



(29) 

where rj s = 4irR^p s /3 is the packing fraction of the 
spheres and Q Sy b is the excess free energy density of a 
pure hard-sphere fluid, a can be considered as the free- 
volume fraction, i.e., the relative amount of the volume 



V that is accessible to the platelets. Besides providing 
access to inhomogeneous density distributions the above 
geometry-based density functional theory for mixtures 
of hard spheres and platelets offers in addition a sys- 
tematic approach to calculate the work W = — ksThia 
required to insert a platelet into a solution of spheres. 
The expression for the free- volume fraction a in Eq. I|29|l 
is equivalent to the result from a recent scaled-particle 
approach |l0j ]. For comparison we note that the corre- 
sponding free-volume fraction for thin rods is given by 
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exp 



kLR 2 Ps 
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(30) 



From the bulk grand potential function all thermody- 
namic quantities can be calculated. Equating the pres- 
sure and the chemical potentials of both species in both 
phases yields the coexisting densities. 

Figure|3(a) displays the calculated phase diagram for a 
binary mixture of spheres and thin platelets for size ratio 
R s / R p — 8/3 as a function of the chemical potential p p of 
the platelets and the number density p s of the spheres. 
The tie-lines are horizontal because of the equality of 
p p of the coexisting phases. The binodal for coexisting 
states is shown, where a sphere-rich and a platelet-poor 
liquid phase coexists with a sphere-poor and a platelet- 
rich liquid phase. The coexistence region is bounded by 
a lower critical point below which only a single stable 
phase is found. For convenience we have introduced the 
dimensionless variable p* = p p — ksT \n(A p / R p ), and 
dropped the star in order to avoid a clumsy notation. 
Figure |3|(b) displays an alternative representation of the 
phase diagram in terms of the number densities of both 
species. The figure illustrates the fractionation of both 
spheres and platelets due to the phase transition. Upon 
increasing the size ratio the critical point shifts to larger 
densities of the spheres. 

The dotted lines in Fig. [31 represent the binodal as cal- 
culated without the term proportional to p 2 in paren- 
thesis on the r.h.s. of Eq. (|29|l which is equivalent to 
considering a binary mixture of spheres and thin rods of 
length L = ttR p + 2R 2 /R S = 1.46 R s as is apparent from 
a comparison of Eqs. (|2T?|l and (|3Tj|) . The phase bound- 
aries and the lower critical point are shifted to smaller 
values of p s . In other words, the number density of the 
spheres in the sphere-poor phase is increased by taking 
into account the last term in parenthesis in Eq. I|29|) . 



IV. BINARY SPHERE-PLATELET MIXTURE 
NEAR A PLANAR HARD WALL 

The density profiles of both components of the binary 
mixture of hard spheres and thin platelets close to a pla- 
nar hard wall are obtained by a numerical minimization 
of the grand potential functional Q with the excess free 
energy functional given by Eq. J2J). We fix the chemical 
potential p p of the platelets and approach the bulk phase 
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FIG. 3: (a) Bulk and surface phase diagrams of binary mix- 
tures of hard spheres of radius R a and thin hard platelets 
of radius R p — 3R s /8 as a function of the chemical poten- 
tial of the platelets p p and the number density of the spheres 
p s - The straight dashed lines are tie-lines illustrating liquid- 
liquid phase coexistence, (b) Phase diagram of the same fluid 
in the density-density plane, where p v is the number density 
of the platelets. In (a) and (b) the solid and open circles de- 
note the bulk critical point and the wetting transition point, 
respectively. Between the wetting transition point and the 
critical point the sphere-rich liquid phase completely wets the 
interface between the hard wall and the sphere-poor liquid 
phase. The dotted curve is the binodal as calculated without 
the last term in parenthesis on the r.h.s. of Eq. 129H which 
is equivalent to the bulk phase diagram of a binary mixture 
of spheres of radius R a and thin rods of length L = 1.46 R s 
[see the main text]. In Fig. 0] density profiles near a hard wall 
are shown along the thermodynamic path indicated by the 
arrow at /i p /(feT) = —4.85 in (a); in (b) this path would run 
parallel to the dashed tie-lines. 



boundary from the sphere-poor side. Upon decreasing p p 
the adsorption behavior changes qualitatively, and it is 
worthwhile to distinguish the following two cases. For 
fjLpKksT) > —4.84 we find that the wall is only partially 
wet by the spheres. The layer thickness of the sphere-rich 
phase forming close to the wall increases continuously, 
but remains finite at coexistence. For /i p /(fcsT) < —4.84 
we observe complete wetting. The transition to complete 
wetting appears to be first order because the excess ad- 
sorptions Ti = Rf dz (pi(z) — pj(oo)) with I = p,r, 



FIG. 4: Equilibrium density profiles of hard spheres 
of radius R s (a) and thin hard platelets of radius 
R p = 3R s /8 (b) in contact with a planar hard wall 
at 2 = as the bulk phase boundary is approached 
along the path indicated by the arrow in Fig. [3] (a). 
The chemical potentials of the spheres are /J s /(feT) = 
25.82862, 25.82866, 25.8287, 25.82874, 25.82877 corresponding 
to p s R 3 a = 0.058057,0.058069,0.058082,0.058094,0.058104 
(from left to right) where the chemical potential at bulk co- 
existence is p s /(k B T) = 25.8287774 so that p s Rl = 0.058107. 
In (a) p 3 (z < R s ) =0 and in (b) p p (z < 0) = 0; the contact 
values at z — R a and z — 0, respectively, are not shown on 
the present scales. 

calculated along the coexistence curve jump to a macro- 
scopic value upon approaching the wetting transition 
point. Figure 0] (a) displays the sphere density profiles 
at ^p/(fcsT) = —4.85 signalling the growth of a thick 
layer of sphere liquid at the wall. The corresponding 
platelet profiles are shown in Fig. 01 (b) and indicate how 
the platelets become more depleted as the sphere-rich 
layer grows. Upon approaching the chemical potential of 
the spheres at bulk coexistence pbs/iksT) = 25.8287774, 
the calculated density profiles at the liquid-liquid inter- 
face become virtually indistinguishable from the ones of 
the free liquid-liquid interface between coexisting bulk 
phases, and the layer thickness diverges logarithmically, 
as expected for the case of complete wetting in systems 
governed by short-ranged forces. With increasing chem- 
ical potential of the platelets, and hence increasing dis- 
tance to the critical point [see Fig. 01(a)], the interface 
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becomes sharper, i.e., it crosses over from one to the 
other limiting bulk value over a shorter distance. The 
wavelength A = 1.72 R s of the oscillations of the density 
profiles close to the wall reflects the size of the spheres. 

Similar to recent studies of wetting in sphere-polymer 
[jj 0| and sphere-rod £| mixtures, we have not been 
able to numerically resolve the prewetting line which 
should emerge tangentially from the coexistence curve 
at the wetting transition. In contrast to those studies 
there are no layering transitions in the partial wetting 
regime for the sphere-platelet mixture. This holds also 
for the aforementioned toy model without the last term 
in Eq. (JIJ which exhibits the same bulk phase diagram 
as the corresponding sphere-rod mixture (see the dotted 
curves in Fig. 0). Taking into account the weight func- 
tion (z, P ) [Eq. (fTp|) ] leads to wetting phenomena of 
sphere-platelet mixtures which are different from those of 
sphere-rod mixtures even if the bulk phase diagrams of 
both systems are identical. Finally, we note that the wet- 
ting behavior discussed above remains unchanged upon 
increasing the size ratio of the platelets and the spheres 
(e.g., Rp = R s /2). 

V. SUMMARY 

We have developed a geometry-based density func- 
tional theory for fluids consisting of hard spheres and 
thin platelets. The bulk and surface phase diagram and 
the density profiles near a planar hard wall are deter- 
mined numerically with the following main results. 

(1) Figure ^ illustrates that from their side view thin 
platelets may be regarded as thin rods and from their top 
view as two-dimensional spheres. On the basis of this 
consideration we have shown that the geometry-based 
density functional theory developed for binary mixture 
of hard spheres and thin rods 0, 0, can be 



consistently extended to the more challenging problem 
of hard spheres mixed with thin hard platelets by intro- 
ducing an additional weight function which characterizes 
the surface of a platelet. The volume accessible to a 
thin platelet of radius R p in the presence of a sphere is 
smaller than the corresponding one of a thin rod of length 
L = 2R P because of the extended surface of the platelet 
(Fig.©. 

(2) The bulk phase diagram exhibits two-phase coexis- 
tence between sphere-rich and sphere-poor phases which 
is bounded by a lower critical point below which a sin- 
gle stable phase is found (Fig.EJ. The phase boundaries 
and the critical point of the corresponding phase diagram 
for a binary mixture of hard spheres and thin rods are 
shifted to smaller values of the density of the spheres 
due to smaller intermolecular interactions between thin 
rods and spheres as compared with those between thin 
platelets and spheres. 

(3) For the mixture near a planar hard wall, a first- 
order wetting transition by the sphere- rich phase occurs. 
In the partial wetting regime no layering transitions are 
found in contrast to recent studies on binary mixtures of 
spheres and non-adsorbing polymers or thin hard rods. 

We have focused on the case of non-interacting 
platelets as regards their mutual interactions, which con- 
stitutes a minimal model for non-spherical particles with 
non-vanishing surface area. With increasing density of 
the platelets, interactions between platelets must be in- 
cluded [22, 123 • The consistent treatment of these non- 
trivial platelet-platelet interactions within a geometry- 
based density functional theory remains as a challenge. 
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